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摘 要 在 实际 运行 过 程 中 ， 槽 式 太阳 能 集 热 器 的 聚 光 过 程 可 能 受到 太阳 光 入 射 角 、 安 装 误差 、 镜 面 误 差 和 跟踪 误差 等 各 种 
实际 因素 的 影响 。 本 文 用 空间 坐标 变换 的 方法 结合 蒙特 卡 罗 光 线 追 踪 法 ， 分 别 计算 出 了 各 种 实际 因素 单独 存在 时 和 多 种 实际 
忆 素 耦合 时 的 槽 式 集 热 器 热流 分 布 。 结 果 显 示 ， 安 装 误 差 和 跟踪 误差 使 聚焦 热流 的 不 均匀 性 增 大 ， 镜 面 误 差 使 聚焦 热流 的 不 
均匀 性 减 小 ， 计 算出 的 实际 因素 耦合 热流 分 布 ， 可 以 为 模式 集 热 器 光学 性 能 和 热 性 能 的 分 析 提 供 依据 。 
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Investigation on coupled optical factors of concentration process in 


parabolic trough collectors 
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Abstract The concentration process of a parabolic trough collector can be affected by incidence angle, installation 
error, reflector slope error and tracking error in operation. The coordinate transformation method and the Monte 
Carlo Ray Trace method are integrated to compute heat flux distributions under individual and coupled optical errors 
respectively in this study. It is shown that installation error and tracking error increase the non-uniformity of the heat 
flux distribution while reflector slope error tends to decrease it. Besides, the heat flux distribution under coupled 
optical errors can provide foundations for thermal and optical analysis of parabolic trough collectors. 
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0 名 | 膏 1.1 热流 分 布 模拟 方法 

模式 太阳 能 集 热 器 (PTC) 是 太阳 能 热 发 电 的 ee 
一 种 重要 方式 ， 在 对 PTC 进行 五 究 时 ， 集 热 管 外 聚 画 出 玻璃 管 。 反射 镜 张口 宽度 、 焦距 和 集 热 器 长 度 
焦 热 流 的 计算 常常 是 第 一 步 。 何 雅 玲 中 使 用 蒙特 卡 分 别 用 w,f 和 L 表示 ， 坐 标 原点 建立 在 吸 热管 的 
罗 光 线 追 踪 法 (MCRT) 计算 了 理想 情况 下 PTC 的 。 进口 截面 圆心 。 了 
聚焦 热流 ， 证 明了 MCRT 方法 的 实用 性 ，MCRT 方 MCRT 方法 的 细节 见 文献 [1]。 在 MCRT 方法 
法 得 到 了 广泛 的 应 用 P23。 与 理想 情况 相 比 ，PTC 的 。 ”中 ， 假 如 入 射 光线 方向 向 量 为 MM ， 反 射 面 的 法 向 


聚 光 过 程 受到 各 种 实际 因素 包括 太阳 光 入 射 角 、 安 量 为 ， 那 么 反射 光线 方向 向 量 P 的 计算 公式 为 : 


装 误差 、 跟 踪 误差 和 镜面 误差 的 影响 。 赵 东明 等 B= 2(N.H)N+H 加 
采用 空间 坐标 变换 的 方式 研究 了 集 热 管 安装 误差 、 i ee | 
跟踪 误差 和 太阳 光 入 射 角 对 PTC 聚焦 热流 的 影响 ， cb els sri 


程 泽 东 等 外 对 镜面 误差 和 跟踪 误差 进行 了 详尽 的 敏 。。 线 或 借 国 法 线 万 向 产生 影响 ， 在 MCRT 的 和 这 中 
感性 分 析 ， 均 发 现 当 这 些 实际 因素 存在 时 ， 珍 焦 热 。 。 过 二 考 问世 个 二 水 。 久 了 研究 各 种 因素 的 影响 只 
流 将 发 生 很 大 的 改变 ， 因 此 需要 对 各 种 实际 因素 存 党 要 改变 集 热 入 的 位 置 ， 或 用 空间 红 标 变换 的 方法 


台风 3 A 3 HZ 一 1 水 uM > 和 7 
在 时 PTC 的 集 热管 外 聚焦 热流 进行 研究 。 计算 新 的 人 或 w 即 可 。 ， 
目前 的 文献 常常 只 计算 出 只 有 一 种 实际 因素 存 吸 热管 
在 时 的 PTC 聚焦 热流 ， 而 实际 应 用 中 , PTC 的 聚 光 RA) 


过 程 可 能 会 同时 受到 多 种 实际 因素 的 影响 。 本 文 先 
计算 各 种 实际 因素 单单 独 存 在 时 PTC 的 聚焦 热流 并 7 
进行 分 析 ， 然 后 在 MCRT 方法 的 基础 上 使 用 空间 坐 
标 转换 ， 计 算出 多 种 实际 因素 耦合 时 的 聚焦 热流 。 


> 
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1 模型 描述 (a) 截面 示意 
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(b) 三 维 视图 

图 1 模式 集 热 器 示意 图 

Fig. 1. Schematic of a PTC (a) cross sectional view (b) 3D 
View 


1. 2 安装 误差 

当 集 热管 的 安装 误差 使 吸 热管 的 位 置 发 生平 移 
时 ， 假 设 管子 圆心 线 在 y 和 = 方向 新 的 坐标 分 别 为 
力 和 二， 则 吸 热管 外 壁 方程 变 为 ; 

(y-ny+(2-a} = i 

”是 吸 热管 的 外 壁 半径 。 集 热管 安装 误差 或 者 
支撑 架 的 形变 也 可 能 会 造成 集 热管 的 位 置 与 焦点 线 
不 平行 。 当 PTC 采用 双 轴 跟踪 系统 时 ， 光 子 被 吸 热 
管 接收 位 置 的 x 坐标 等 于 镜面 上 入 射 点 的 x 坐标 
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(用 xo 表示 )， 所 以 只 要 知道 吸 热管 中 心 线 的 坐标 
方程 ， 光 子 到 达 吸 热管 的 位 置 仍 然 可 以 求 出 来 。 

简单 起 见 ， 假 设 吸 热管 的 进口 截面 圆心 坐标 为 
(0,0,0)， 出 口 截面 圆心 坐标 为 (L,yj,z7)， 那 么 在 反射 
光线 所 在 的 垂直 于 焦点 线 的 平面 内 ， 吸 热管 外 壁 的 
坐标 方程 为 : 


2 过 
[| + | = 六 
G3) 


将 吸 热管 新 的 坐标 方程 (2) 或 (3) 代 入 MCRT 算 
法 中 ， 保 持 其 他 部 分 不 变 ， 即 可 求 得 有 集 热 管 安装 
误差 的 管 外 不 均匀 热流 分 布 。 

对 于 单 轴 跟 踪 系 统 ， 由 于 入 射 角 的 存在 ， 反 射 
光线 和 集 热 管 外 壁 的 交点 不 易 确定 ， 情 况 更 为 复 
杂 ， 本 文 不 作 研 究 。 

1.3 太阳 光 入 射 角 

PTC 对 太阳 位 置 的 跟踪 有 单 轴 跟踪 和 双 轴 跟踪 
两 种 ， 只 有 双 轴 跟踪 系统 才能 使 太阳 光线 保持 直射 
入 射 ， 但 是 由 于 双 轴 跟踪 系统 结构 复杂 ， 工 程 实际 
中 PTC 一 般 采 取 单 轴 跟 踪 系 统 ， 这 种 情况 下 入 射 太 
阳 入 射 角 往 往 不 为 0"， 如 图 2 所 示 。 

太阳 光线 4B 以 入 射 角 0 照射 到 反射 镜 边缘 上 
的 点 B， 然 后 被 反射 到 吸 热 管 上 的 点 C， 距 离 CD 
即 为 经 过 B 点 的 太阳 光线 不 能 照射 到 的 长 度 。 文 献 
上 求 出 CD 的 最 大 值 和 最 小 值 分 别 为 


Ww 
钱 、… 二 +——— ltan0 
max [ 号 | 


(4) 
Loin > ftan0 


(5) 
Lmax 是 PTC 末端 损失 效应 能 够 影响 到 的 最 大 距 
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离 ， 即 与 末端 DD 距离 大 于 Zax 的 区 域 总 的 聚焦 热流 
不 变 ， 而 与 末端 DD 距离 小 于 Lwin 的 管 段 则 不 能 接收 
到 反射 后 的 太阳 光 。 


六 太阳 


图 2 太阳 光 斜 射 入 射 图 
Fig. 2. Diagram of an incident ray on the Sun-facing end ofthe 
reflector 


图 3(a) 中 光线 1 是 MCRT 法 计算 出 来 的 直射 入 
射 光 线 ， 光 线 2 是 入 射 角 为 0 的 和 斜 射 入 射 光 线 。 图 
3(b) 中 的 坐标 系 O-xy2 是 坐标 系 O-xyz 绕 y 轴 旋 转 
角度 9 后 得 到 的 ， 光 线 1 在 此 坐标 系 中 的 坐标 就 是 
光线 2 在 坐标 系 O-xyz 中 的 坐标 ， 用 17 在 表示 ， 
则 ; 


M'=A:M (6) 
其 中 4 是 从 坐标 系 O-xyz 到 坐标 系 O-xy2 的 空 
间 坐 标 变换 矩阵 ; 


coSO_ 0 SmnO 
4 = 0 1 0 
-SInO 0 cosO 


(7) 
将 得 出 的 代入 公式 (1) 中 ， 保 持 其 他 部 分 不 
变 ， 即 可 求 得 太阳 光 入 射 角 不 等 于 0° 时 的 热流 分 


(b) 空间 坐标 变换 


3 单 轴 系 统 太 阳光 入 射 角 坐 标 计算 
Fig. 3. The calculation of the incident direction of a single-axis 
PTC (a) different incident conditions where ray 1 is the vertical 
incident ray and ray 2 is the oblique incident ray (b) the 
coordinate transformation for the incidence angle 0 


1.4 跟踪 误差 
跟踪 误差 是 由 于 跟踪 不 准确 导致 的 集 热 器 位 置 


PE 


的 偏差 ， 如 图 


4(a) 所 示 。 由 于 跟踪 误差 的 存在 ， 集 


热 器 1 相对 于 正确 的 位 置 〈 集 热 器 2) 逆 时 针 旋 转 


了 8 角度 。 如 图 4(b) 所 示 ， 空 间 坐标 变换 矩阵 等 


0 是 标准 方差 ， 也 是 镜面 误差 角度 。 将 得 出 的 
WW 代入 公式 (1) 中 ， 保 持 其 他 部 分 不 变 ， 即 可 求 得 
有 镜面 误差 的 热流 分 布 。 
1.6 实际 因素 耦合 

对 于 安装 误差 ， 与 理想 模型 相 比 ， 吸 热管 坐标 
发 生 了 变化 ， 用 式 (2) 或 式 (3) 来 考虑 其 影响 即 可 。 对 
于 跟踪 误差 、 镜 面 误差 、 太 阳光 入 射 角 等 实际 因 


于 
1 0 0 
有 =|0 cosp sinf 
[10 -sn cosp (8) 
则 修正 后 的 入 射 光线 坐标 27 为: 
M’=B:.M (9) 


将 得 出 的 代入 公式 (D) 中 ， 保 持 其 他 部 分 不 
变 ， 即 可 求 得 


素 ， 可 以 用 以 下 公式 计算 反射 光线 方向 ， 
P=-2(N .MT ) N+ MM 


(13) 


其 中 WV’ 和 MW 分 别 是 考虑 了 跟踪 误差 、 镜 面 误 


差 和 太阳 光 入 射 角 等 因素 的 实际 反射 镜 法 


线 和 实际 


跟踪 误差 的 热流 分 布 。 
这 


M (CO) 


入 射 光 线 坐标 向 量 。 


需要 注意 的 是 ， 由 于 矩阵 乘法 不 满足 交换 律 ， 


所 以 坐标 变换 的 顺序 不 是 随意 的 。 


例如 ， 当 同时 考 


虑 入 射 角 影 响 和 跟踪 误差 时 ， 实 际 入 射 光 线 坐 标 11; 
为 : 


M'=B.(4:M) 


(14) 


(b) 空间 坐标 变换 
图 4 跟踪 误差 及 其 坐标 变换 
Fig. 4 Tracking error and its coordinate transformation (a) 
diagram of tracking error (b) the coordinate transformation for 
the tracking error f 


1.5 镜面 误差 
镜面 误差 考虑 的 是 反射 面 微 观 的 缺陷 或 宏观 的 


采用 MM?=B.(4: 必 ) 或 VI=4.(8. 放 得 到 的 
者 果 并 不 一 样 。 算 阵 4 把 直射 光线 转变 为 斜 射 入 射 
光线 ， 和 矩阵 B 考虑 跟踪 误差 ， 考 上 处 到 实际 过 程 ， 应 
该 采取 2M'=B.(4. 斩 ) 来 计算 跟踪 误差 存在 时 单 
跟踪 系统 的 热流 分 布 。 
2 结果 与 讨论 
下 面 利用 空间 坐标 变换 和 矩阵， 对 安装 误差 、 太 
阳光 入 射 角度 等 实际 因素 分 别 进行 模拟 和 分 析 ， 最 
后 将 多 种 实际 因素 耦合 在 一 起 计算 。 因 为 管 外 热流 
分 布 的 相对 大 小 与 直射 辐射 强度 没有 关系 ， 所 以 分 
析 结 果 时 管 外 热流 取 的 是 当地 聚焦 比 (LCR)。 本 文 
模拟 时 采用 LS-3 型 模式 太阳 能 集 热 器 作为 物理 模 
型 ， 几 何 尺 寸 等 参数 见 表 1。 关 于 其 他 型 号 集 热 器 

的 结论 可 以 通过 用 相同 的 方法 分 析 得 到 。 

表 1 LS-3 模式 集 热 器 数据 

Table1 Parameters of the LS-3 PTC 


NS 


变形 等 原因 所 造成 的 反射 误差 。 镜 面 误差 包含 的 种 - 夫 攻 数值 
类 很 多 , 把 所 有 的 镜面 误差 着 加 起 来 后 可 以 用 实际 吸 热管 半径 /m 0.035 
的 与 理想 的 镜面 法 线 方向 的 偏差 来 代 蔡 m， 认 为 此 焦距 血 1 
偏差 角度 在 x 方向 和 方向 分 解 后 得 到 的 两 个 分 量 。 ”反射 镑 张口 宽度 加 576 
分 别 符合 相同 的 高 斯 分 布 四 ， 镜 面 实际 法 线 坐 标 向 。 玻璃 管 透射 率 095 
量 浆 可 以 由 理想 的 镜面 法 线 方向 不 经 过 两 次 坐标 吸 热管 吸收 率 0.96 
变换 得 到 : 反射 镜 反 射 率 0.93 
Saad (10) ”2 1 安装 误差 的 影响 
日 八 = 过 、 AS 《4 。 E34 
PO a 方向 和 2 方向 安装 误差 的 定义 分 别 为 中: 
人 人 ， Y=yyw 和 2Z=zyf。 图 5 中 显示 的 是 平移 安装 误差 下 
Ge oo | 的 热流 分 布 。 集 热管 安装 误差 达到 毫米 数量 级 时 ， 
Varo \20 (11) 管 外 热流 形状 就 已 发 生 明 显 改变 。 图 6 和 图 7 比较 
了 两 方向 安装 误差 都 存在 时 y 或 z 方向 安装 误差 对 


管 外 热流 分 布 的 影响 。 图 6 保持 Z 不 变 改变 7 的 大 


小 ， 热 流 分 布 形状 的 改变 非常 大 ;图 7 保持 了 不 变 
改变 Z 的 大 小 ， 热 流 分 布 形状 的 改变 相对 较 小 。 因 


此 ,与 z 方 向 安装 误差 相 比 ，y 方向 安装 误差 的 影 


响 大 。 所 以 ， 在 集 热管 的 安装 过 程 中 要 格外 注意 减 


少 在 y 方 向 的 安装 误差 。 
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(b)y 方向 安装 误差 


图 5 


平行 安装 误差 下 的 热流 分 布 


Fig. $5 Flux distributions for installation errors (a) z direction 
installation errors (b) y direction installation errors 
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(b) 2=0.5% 
图 6Z 相 同时 不 同 了 对 聚焦 热流 的 影响 
Fig. 6 Flux distributions for different Y when Z is constant (a) 
Z=0.2% (b) 2=0.5% 
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(b) 二 0.5% 
7 了 相同 时 不 同 Z 的 影响 
Fig. 7 Flux distributions for different Z when 了 is constant (a) 
二 0.2% (b) 二 0.5% 


图 8 显示 的 是 不 平行 安装 对 热流 分 布 的 影响 。 
集 热 管 进口 截面 圆心 坐标 为 (0,0,0)， 出 口 截 面 圆心 
坐标 为 (4;y1,z/)， 其 中 yj 和 zj 分 别 对 应 天 0.1% 和 
Z=0.1%。 虽 然 图 中 热流 分 布 在 轴线 上 是 不 均匀 的 ， 
但 可 以 看 成 无 数 个 平移 安装 误差 的 热流 分 布 连 续 排 
布 而 成 。 


| 


图 8 出 口 截 务 


| 天 0.1%，Z-0.1% 热 流 分 布 
Fig. 8 Flux distributions in case of 7=0.1% and 2=0.1% for 


outlet of the absorber tube 


2.2 入 射 角 的 影响 


图 9(a) 显 示 的 是 直射 辐射 强度 相同 时 双 轴 
和 单 轴 跟踪 系统 热流 


绩 


沿 管 长 的 分 布 ， 纵 轴 


收 的 总 能 量 。 对 于 
Lmax 和 Lmin 的 位 置 


受 末端 损失 效应 影响 区 


的 总 的 聚焦 热流 等 


认为 是 不 变 的 。 从 区 


端 ， 吸 热管 接收 到 凋 


跟踪 的 系统 ， 图 中 画 出 了 


可 以 看 至 


太阳 光 以 相等 的 直射 辐 
直射 入 射 时 接收 的 聚焦 热流 ， 且 在 轴 向 上 仍然 可 以 


为 截面 吸 


| MCRT 方法 模拟 的 结 
MR 
对 太阳 光 斜 射 入 射 时 热流 分 布 模拟 的 准确 性 。 


区 域 ， 集 热管 每 个 截面 接收 


虽 射 强度 


网 有 和 
i ln 


大 ， 最 终 成 为 周 向 不 均匀 、 轴 向 上 不 变 的 聚 


流 。 
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图 9 太阳 光 30°* 入 射 角 热流 分 布 


Fig. 9 Flux distributions at 30° incidence angle (a) axial 


distribution (b) 3D distribution 


CR 

T 60.00 
50.00 
40.00 
30.00 
20.00 
10.00 
0.000 


2.3 镜面 误差 和 跟踪 误差 影响 

图 10 显示 的 是 镜面 误差 对 聚 光 过 程 的 影响 。 
从 图 10(a) 可 以 看 出 ， 周 向 热流 分 布 仍然 是 对 称 的 ， 
镜面 误差 的 作用 是 使 聚焦 热流 由 不 均匀 趋向 均匀 ， 
这 种 作用 有 利于 减 小 管 壁 温差 ， 且 镜面 误差 越 大 越 
明显 。 图 10(b) 显 示 的 是 镜面 误差 角度 对 未 受 末端 损 
失 影响 区 域 的 拦截 系数 的 影响 ， 当 镜面 误差 大 于 
0.12° 时 ， PTC 的 拦截 系数 小 于 1， 应 该 避免 这 种 
情况 。 
图 11 显示 的 是 跟踪 误差 对 聚 光 过 程 的 影响 。 
从 图 11(a) 可 以 看 出 ， 跟 踪 误 差 使 太阳 光 偏 离 垂直 于 
反射 镜 张口 平面 的 方向 入 射 ， 增 大 了 热流 的 不 均匀 
性 。 图 11(0) 显 示 的 是 跟踪 误差 角度 对 未 受 末 端 损失 
影响 区 域 的 拦截 系数 的 影响 ， 当 跟踪 误差 大 于 0.6° 
时 ， 拦 截 系数 小 于 1。 


60 


50 上 


40 上 


0 pe Re 
0 60 120 180 240 300 360 
口 ] 


(a) 周 向 热流 分 布 


0.9 


0.8 


0.7 . l 
0.00 0.15 0.30 0.45 


] [ 
(b) 镜面 误差 对 拦截 系数 的 影响 
图 10 镜面 误差 对 聚 光 过 程 的 影响 
Fig. 10 Effects of reflector slope error on concentration process 

(a) circumferential heat flux distributions (b) variation of the 
intercept factor 
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Fig. 11 Effects of tracking error on concentration process (a) 
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图 12 太阳 光 入 射 角 30*， 玉 0.1%，2Z=0.1%， 跟 踪 误 差 
0.2*， 镜 面 误差 0.1" 时 的 热流 分 布 
Fig. 12 Flux distributions in case of 30° incidence angle， 


=0.1%, Z=0.1%, 0.2° tracking error and 0.1° reflector Slope 


error (a) 3D distribution (b) circumferential distribution 


图 12 显示 了 多 种 实际 因素 耦合 时 的 热流 分 
布 ， 太 阳光 入 射 角 为 30°， 
误差 为 0.2"， 镜 面 误 差 为 0.1°。 


差 角 度 的 大 小 不 一 ， 所 以 管 外 聚焦 热流 的 形状 变化 
多 样 ， 可 以 计算 出 集 热 管 外 实际 因素 看 全 合 的 热流 分 
布 ， 从 而 为 管内 三 维 温度 场 的 计算 提供 依据 。 


下 0.1%， 二 =0.1%， 跟 踪 
图 12(b) 显 示 的 是 未 
受 末端 损失 影响 区 域 的 周 向 热流 。 由 于 各 种 实际 误 


3 经 
角 、 


论 
本 文 分 别 分 析 了 集 然 管 安装 误差 、 太 阳光 入 射 
镜面 误差 和 跟踪 误差 等 实际 因素 对 PTC 热流 分 布 


的 影响 ， 计 算出 了 多 种 实际 因素 耦合 时 的 PTC 热流 分 


布 ， 


得 出 了 以 下 结论 ; 
1) 集 热管 安装 误差 达到 毫米 数量 级 时 ， 育 焦 热 


流 


状 就 已 发 生 明 显 改 变 ，? 方向 的 安装 误差 比 z 


方向 对 热 
外 注意 减少 在 y 方 向 的 安装 误差 ; 


\ 流 分 布 影响 更 大 ， 集 热管 安装 过 程 中 要 格 


2) 在 末端 损失 影响 区 域 ， 吸 热管 接收 到 的 聚 


焦 能 量 沿 轴 向 从 零 开 始 连 续 地 增 大 ， 在 此 区 域 以 


集 热管 每 个 截面 总 的 聚焦 热流 等 于 太阳 光 以 相 


等 的 直射 辐射 强度 直射 入 射 时 总 的 聚焦 热流 ， 且 在 


轴 向 上 可 以 认为 不 变 ; 


3) 镜面 误差 的 作用 是 使 聚焦 热流 由 不 均匀 趋 


向 均匀 ， 这 种 作用 有 利于 减 小 管 
越 大 越 明 显 ， 而 跟踪 
ee 


跟踪 误差 大 于 0.6" 时 ，PTC 的 拦截 系数 小 于 1; 


ee 
热流 的 不 均匀 


吴 差 增 大 至 


(4) 由 于 各 种 实际 误差 角度 的 大 小 不 一 


管 外 实际 因素 耦合 的 聚 
热 性 能 


[5] 


[7] 


[8] 


可 以 计算 出 集 热 
聚焦 热流 ， 从 而 为 光学 性 
的 分 析 提 供 依 据 。 
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